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The repulsive Lieb-Liniger model can be obtained as the non-relativistic limit of the Sinh-Gordon 
model: all physical quantities of the latter model (S-matrix, Lagrangian and operators) can be put in 
correspondence with those of the former. We use this mapping, together with the Thermodynamical 
Bethe Ansatz equations and the exact form factors of the Sinh-Gordon model, to set up a compact 
and general formalism for computing the expectation values of the Lieb-Liniger model both at zero 
and finite temperatures. The computation of one-point correlators is thoroughly detailed and when 
possible compared with known results in the literature. 

PACS numbers: 67.85.d, 05.30.Jp, 02.30.1k, 03.75.Hh 

I. INTRODUCTION 

The physics of one-dimensional interacting bosons is well captured by the Lieb-Liniger (LL) model [TJ. 
Despite its deceptive simplicity, this model has become a paradigmatic example of quantum integrable 
systems since it proved to have a remarkable richness [2]: its Bethe Ansatz equations, for instance, can 
be explicitly derived and used to study its equilibrium properties at zero and finite temperatures [3] . The 
explicit analysis of the weak to strong coupling crossover of this model has also set a precise benchmark 
for approximate many-body techniques [JJ [2] . The efforts done for the computation of the correlation 
functions of the LL model have also greatly stimulated the development of new and general formalisms, 
such as the quantum inverse scattering method [1] or the bosonization approach [5]. 

Nowadays a renewed interest in the LL model has been triggered by its accurate experimental realization 
[SHn]: in quasi one-dimensional traps, the excitations in the transverse directions are effectively frozen 
and, moreover, the coupling of the ultracold bosons to the external environment can be made very weak 
|10) . These recent experimental advances have opened new perspectives in the field of strongly correlated 
quantum systems: in such a highly controllable set-up it is, in fact, possible to thoroughly investigate 
problems of general nature concerning quantum extended systems, such as the dynamics of integrable 
systems in the presence of small non-integrable perturbations (e.g., three-body interactions and/or a 
weak external trapping potential), the issue of thermalization in quantum integrable and non-integrable 
systems [12] and the behavior of various susceptibilities and response functions. 

The key quantities to answer all these questions are the correlation functions of the LL model. Despite 
the integrability of the model, their explicit computation turned out to be an interesting theoretical 
challenge . For this reason many different approaches have been developed over the years to tackle 
different aspects of this difficult problem: a partial list includes bosonization (which gives the correct long- 
distance behavior of correlation functions) [5] [T3] [14] , quantum Monte Carlo simulations [15] , algebraic 
Bethe Ansatz T6], analytical- numerical methods based on the exact Bethe Ansatz solution [TjJ ITS] . 
Bogoliubov weak- [19] and strong-coupling methods [20] . renormalization group |21] , numerical results 
using stochastic wave-functions [22] and imaginary time simulations [231 [21] . Exact results based on the 
Yang- Yang equations and the Hellmann-Feynman theorem arc presently available for local two-body 
correlations 20,25, while the local three-body correlations were determined at zero temperature in [26] . 

A new method was proposed recently to compute expectation values in the LL model [27J : it exploits 
a different route from all the previous approaches for it is based on an exact mapping between the 
non-relativistic LL model and the relativistic integrable Sinh-Gordon (sh-G) model. This proposal not 
only provides a remarkable simplification of the problem but applies equally well both at zero and finite 
temperatures. In a nutshell, the logical steps on which the method is based are the following: 

1. Due to the relativistic invariance and quantum integrability of the Sh-G model, it is possible to set 
up functional equations [211 I2H] for the matrix elements of its local operators on the asymptotic 
states - known as form factors - and to find their exact solutions |30[ 131] . 

2. The finite temperature and finite density effects of the Sh-G model can also be controlled by solving 
the Thermodynamical Bethe Ansatz (TBA) equations [321 133] • 
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3. In a proper non-relativistic limit of the sh-G model, all quantities of this theory - S'-matrix, La- 
grangian, form factors, TBA equations and so on - reduce to those of the LL model and therefore 
can be used to establish an explicit mapping between the two models. In particular, from the exact 
and known expressions of the form factors of the sh-G model we can explicitly obtain the matrix 
elements of the operators of the LL model we are interested in. 

4. To actually compute the LL correlation functions we have to take into account another aspect of 
the problem, that is, that the LL correlation functions refer to the ground state of the gas at a 
finite density and at a finite temperature, while those of the sh-G model refer to the vacuum state 
(i.e. the state without any particles). This apparent difficulty can be, however, readily overcome by 
using the LeClair-Mussardo formalism [34] that, as a matter of fact, is based on the same quantities 
mentioned above, i.e., the form factors and the TBA equations. 

In this paper we provide the details of the mapping between the sh-G and the LL models first established 
in [27], presenting additional results and focusing values. The paper is organized as follows: in Section 
|LT| we introduce the LL and sh-G models and we recall their TBA equations. In Section III we discuss 
the form factors of the sh-G model and their correspondence with the operators of the theory. In Section 



IV we explore in detail the non-relativistic regime of the sh-G model and we explain how the LL model 
emerges in this limit: the mapping between the two models is discussed at the level of their S'-matrices, 
Lagrangians and TBA equations. In Section [V] we describe the method of calculating local correlation 



functions in the LL model, providing a detailed derivation of the main result, formula (100). Then we 
present the explicit computation of LL one-point correlators at zero and finite temperatures, together 
with their comparison with known exact and approximate results available in literature. Our conclusion 
and outlooks are given in Section |VT) Supplementary material is presented in the Appendices, where we 
also list the non-relativistic limit of the sh-G form factors for the first few cases. 



II. THE MODELS 

In this section we recall the main properties of the LL and sh-G models and we also discuss the expressions 
of their free energy and T — ground-state energy obtained by the TBA equations. 

A. Lieb— Liniger Hamiltonian 

The LL Hamiltonian describes N non-relativistic bosons of mass m in one dimension, interacting via a 
two-body repulsive (5-potential: 

fc2 N o2 

i—l L i<j 

For cold atomic gases the quantity A > in the Hamiltonian can be determined in terms of the pa- 
rameters of the three-dimensional Bose gas in the quasi-one-dimensional limit [10 . The effective coupling 
constant of the LL model is given by the dimensionless quantity 

2mA 

7 = ~T2 — ) 2 

where n = N/L is the density of the gas (L is the length of the system). The limit 7 <C 1 is the weak 
coupling limit and in this regime it is known that the Bogoliubov approximation - obtained by linearizing 
the Gross-Pitaevskii equation - gives a good estimate of the ground-state energy of the system pQ . For 
large 7 one approaches the Tonks-Girardeau limit [35] and recently it has become an interesting question, 
also from the experimental point of view, to study the crossover between the two regimes [2J. In the LL 
model temperatures are usually expressed in units of the quantum degeneracy temperature 

h 2 n 2 

kbTd = — — , (3) 
2m 



in the following we use the scaled temperature r = T/Td. 
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In second quantized formalism the non-relativistic field theory describing bosons interacting via a 
^-potential is defined by the Hamiltonian [4] 

where the complex Bose field ip(x,t) satisfies the canonical commutation relations 

[tp(x,t),^(x',t)] = 5(x-x') , [tl)(x,t),tp(x',t)] =0. (5) 
The Lagrangian density associated to the field theory Hamiltonian T-L is 

.h (.*ty dtf A » 2 s^t^ wt/t// 

Restricting to the subspace of the Hilbert space where the number of particles N is fixed, the equation 
of motion for the field ip translates into the eigenvalue equation of the LL many-body Hamiltonian ([I]) : 

Hxn{xi,---,x n ) = E n xn{xi,...,x n ). (7) 

As shown by Lieb and Liniger in their original paper pQ, the eigenvalue problem Q can be solved in 
terms of a Bethe Ansatz. Let us recall the main steps because they will lead us to the definition of the 
S'-matrix of the LL model, a basic quantity of its dynamics and of our following discussion. One can easily 
see that Eq. ^ is just a free Schrodinger equation in the domain where the coordinates of the particles 
are all distinct. However, we have to enforce the usual boundary conditions for a (5-function potential, 
that is, the discontinuity of the gradient of the wave function when two coordinates coincide 

L " d^) XN ^=^ = AXivU^ . (8) 

If we denote by i?i the subset of the configuration space where X\ < X2 < • • • < xjq , the solution of the 
equations in R\ is given by the Bethe wave function 

X n(xi,x 2 , ...,x n ) = J2 <P) ^ Ef=I P(ki)Xj > (9) 
p 

where ^ p denotes a sum over permutations of the momenta {fci, . . . ,k n } that characterize the state. 
For configurations outside R\ the solution is easily obtained using the symmetry of xn with respect to 
the Xi. The coefficients in the sum are related by the boundary conditions Q. For the permutations 
P : (k, I, k a3 , . . . , k aN ) and Q : (l,k, k a3 , . . . , k aN ) the relation between the corresponding coefficients is 

h — l — j2m \ 

a(Q) = ^— a(P) . (10) 

k-l + i^X ( ' 

Hence the wave function gets multiplied by the factor a(Q)/a(P) whenever two particles with momenta 
Pi = k and p 2 = I are exchanged. This exchange is equivalent to a scattering process of the two particles 
and therefore the two-body 5-matrix of the Lieb-Liniger model is expressed by 



■n 4 2 m \ 



S LL ( P ,A)= ■ . * (11) 



where p = p\ — p2 is the momentum difference. 

Given the integrability of the model, all its physical properties can be essentially derived from its two- 



body S'-matrix ( 11 ). This quantity, for instance, determines the thermodynamics of the model, as shown 
originally by Yang and Yang [3]. In the limit N — > oo, L — > oo with the density n fixed, the discrete 
energy levels of the system get encoded in an energy level density function p(p) and in the density p^ T '(p) 
of the occupied levels. Notice that we are going to use a tilde"" for the quantities in the non-relativistic LL 
TBA, while the corresponding quantities in the TBA for the sh-G model will be later denoted without 
this tilde. The ratio between the two densities p and p^ defines the pseudo-energy i(p) through the 
relation 



(12) 
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and this quantity, together with the densities, satisfies the coupled set of integral equations 

oo 

2tt p(p) = l+ ! dp'$(p-p')pW(p'), (13a) 



zip) 



jr 



fi 



kviT 2rakvT 



2^ 



^(p-j/)log(l + e-^) , 



(13b) 



(13c) 



where p is the chemical potential, T is the temperature and fee is the Boltzmann constant. The kernel 
(pip) that drives all the integral equations entirely follows from the 5*-matrix (11) 



(pip) = -i—\ogS LL {p) = AhrnX 

dp h 2 p 2 + 4m 2 A 2 



(14) 



Once the TBA integral equations (13) are solved, the ground state energy E and the free energy F of 
the system are expressed as 



E 
L 



dp^ip). 



F . k B T 



dp log ( 



1 + e 



-s(p) 



(15a) 
(15b) 



At zero temperature the energy level density gets a compact support, i.e. it is different from zero only 
on an interval (which we denote by [—B,B]) and, correspondingly, the TBA equations simplify as 

B 

2npV (p) = i + / dp' (pip- p') p« (p') , (16a) 

B 

ioip) = -V+7^+ f %viP~p') iotf) , (16b) 



-B 



where ioip) = limr->o k^T eip) and the boundary value B is determined by the normalization condition 

B 

n = f p (r) (p)dp. (17) 
—B 

The TBA equations of the sh-G model will be described in the next section where the physical meaning 
of the pseudo-energy will also be discussed. In Section |IV C| we show that the TBA equations and the 
excitation spectrum of the LL model can be obtained in a proper limit from those of the sh-G model. 



B. The Sinh— Gordon model 



The sh-G model is an integrable relativistically invariant field theory in 1 + 1 dimensions defined by the 
Lagrangian density 
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where 4> — (f)(x, t) is a real scalar field, too is a mass scale and c is the speed of light. The parameter mo 
is related to the physical (renormalized) mass M of the particle by [36] 

m% = M 2 ^- (19) 
sm(7raj 

The explicit presence of the speed of light c will help us in studying later the non-relativistic limit of this 
theory (see Section IV). Despite the relativistic nature of the sh-G model, its integrability (supported 
by the existence of an infinite number of conservation laws) implies the absence of particle production 
processes and that its n-particle scattering amplitudes are purely elastic. Moreover, they factorize into 
n(n — l)/2 two-body S-matrices. The energy E and the momentum P of a particle can be written as 
E = Mc 2 cosh^, P = Mc sinh where is the rapidity. In terms of the particle rapidities, the two-body 
S'-matrix is given by [37]: 

sinh0- i sin(aTr) 

<Ssh-G(0, a) = x , (20) 

siiuiS + I sm(a7r) 

where 9 is the rapidity difference and a is the dimensionless renormalized coupling constant 

8tt + he g 2 K ' 

As in the LL model, the two-body S'-matrix of the sh-G model fully encodes its physical properties, 
in particular its thermodynamics. Its derivation is quite similar to the one of the LL model, the only 
difference being the relativistic kinematics [35J [33J . Let us briefly discuss the TBA equations of the sh-G 
model both at finite temperature and at finite particle density n. Note initially that, although the sh-G 
model is a relativistic theory, its quantum integrability implies the conservation of the number of particles 
and therefore it makes sense to associate a wave-function to the iV-particle state. Therefore, the starting 
point of the TBA approach is the quantization of the rapidities of the TV-particle state on an interval L 
with periodic boundary conditions, given by 

JV 

McL sinh 9 t +H^2 X(9i ~ 6j) = 2ttM^ , (22) 

where the Mi are (positive or negative) integers and = — * log <Ssh-G(#) is the phase shift of the two- 
body scattering process of the sh-G model. In the thermodynamic limit (N — > oo, L — > oo, N/L = n = 
fixed) we introduce for the left hand side of Eqs. (22) the quantity 

oo 

J (9) =Mc sinh 6 + 2%K [ — x (9 - 9')p {r) (9') , (23) 

J 2-7T 

— oo 

and by differentiating it equations (22) turn into the integral equation 

oo 

P{&) = S C0Sh6 ' + J f%(0-0')P W (0')- (24) 

— OO 

In the equations above 

is the density of states, p^ T '(9) is the density of the occupied states (both per unit length) and <p(9) is the 
derivative of the phase shift 

i p(0) = ^p- = -i^ogS sh . G (9), (26) 



The pseudo-energy e(9) is introduced as in ([12]) 

= l + e eW . (27) 



P M(0) 
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By minimizing the free energy 
F\p,pV]=E-TS = 

L [ M cosh 6 p (r) (0) d6-LT I p logp - p M logp (r) - (p - p (r) ) log(p - p (r) ) d9 (28) 



with respect to the densities p(9) and p^(6), with the constraint (24), one arrives at the TBA equation 
for the pseudo-energy 
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(29) 



Once this integral equation for e{&) has been solved, the densities p(9) and p^ r \9) are extracted from 
equations ( 24 1 and ( 27 ) while the chemical potential is fixed by the constraint 



= / d0p«(0) 



(30) 



Any other thermodynamic quantity can be calculated from the free energy (per unit length) of the system 
determined from the above minimum principle and finally expressed by the formula 



f= F = k B T 
1 L 2nh 



d9Mccosh9 log (1 + e 



pn . 



(31) 



Let us consider now the T — > limit of the TBA equations. Due to the non-zero chemical potential, we 
can assume that the function e(9) changes sign at the rapidity values —9* and 9* so that it is negative 
on the interval I = (— 9*, 9*), zero at the limiting points and positive everywhere else. For T — > 0, e{9) 
becomes largely negative on the interval I and largely positive outside. Then from Eq. ( [27] ) we see that 
on the interval (—9*, 9*) we have p( T ' = p (and a filling fraction equal to 1), while outside p^ = (and 
a filling fraction equal to 0). The boundary value 9* is determined by the condition 



n = / p^ (jp) dp . 



(32) 



For -9* < 9 < 9* equations (24) and (29) become 



pM(0) = ^coshfl 
zirn 



d6' 
2^ 



p{9-9')p^\9'). 



(33a) 



e o (0) = Mc 2 cosh 9 - p 



d0' 



V (9-9')e {9'), 



(33b) 



where e${9) — limy-i-o kBTe(9). 

Let us pause here to comment on the physical meaning of e(9). If we raise one of the numbers Afi in 
the Bethe equation that corresponds to a rapidity 9 i: to a larger value Af!, which will correspond to a 
rapidity 9', then the change in the energy of the system is 

AE = k B T(e(9 / ) - e(9)) . (34) 

This shows that e(9) describes the energy of the excitations over the ground state, and it gives the dressed 
energy of the quasi-particles (hence the name pseudo-energy) . If we demand that the ratio of the densities 
take the usual form 

gW(g) _ 1 
p(9) ~ e (-E(e)-M)/fc B T + i 

then the excitation energy is fixed to be 

E{6) = k B T e{9) + p . (36) 



(35) 
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III. FORM FACTOR EXPANSION FOR ONE-POINT CORRELATORS 



Expectation values in an integrable relativistic field theory are conveniently expressed in terms of the 
so-called form factors. 

The form factors we are going to use in this article are the ones related to a quantum field theory, which 
in principle are different from the quantities that share the same name in the Bethe ansatz solution of 
integrable models. In the latter context they are defined as matrix elements of operators between exact 
Bethe ansatz states, while here we use instead the basis in the Hilbert space consisting of multiparticle 
scattering states. However, these quantities turn out to be closely related [38 . In this section, for the sake 
of completeness, we first review the definition and the main properties of the form factors in a general 
relativistic field theory. Since our final goal is to treat the sh-G model which is a theory with a single 
type of gapped excitations and multi-particle states without any bound states, we focus our attention on 
theories of this type. (A well-known example of a theory with different types of excitations and bound 
states is the Sine-Gordon model [31].) Later we present the explicit expressions of the form factors of the 
sh-G model. 



A. Basic properties of form factors 

Consider a local operator 0(x, t). Using the translation operator U = e~ tp ' iX , where t — x° and x = x 1 , 
we can always shift this operator to the origin 0(x,t) = WO(0,0)U. The matrix elements of 0(0,0) 
between the vacuum and a set of n-particle asymptotic states are called the n-particle form factors of 
this operator (see Fig. [T]) 

F°(9 1 ,9 2 ,...,9 n ) = (O|0(O,Q)|0i,...A). (37) 

For an operator of spin s, relativistic invariance implies that under a simultaneous shift in the rapidities 
its form factors behave as 

(0x + A, ...,#„ + A) = e sA F^(9 1) . . . , 9 n ) . (38) 

This equation indicates that the form factors of a scalar operator depend only on the differences of 
rapidities, 6ij = 9i — 9j. A generic matrix element of the operator 0(0,0) can be expressed in terms of 
its form factors by using the translation operator and the crossing symmetry, which is implemented by 
an analytic continuation in the rapidity variables 

(6 1 , . . . ,e n \O(0, 0)101, . . . ,p m ) = F% +m (p x , ...,[3 m ,6 x -m,...,e n -m). (39) 

(If Pi = 6j for some i and j, this formula gets modified by contact terms.) Hence, the knowledge of all 
form factors of an operator is equivalent to the knowledge of the operator itself, (i.e., how it acts on any 
state of the theory) . 

The form factors satisfy a set of functional and recursive equations, which for integrable models makes 
it possible to find in many cases their explicit expressions (for a review, see [28] [40] ) . For a scalar operator 
the functional equations (known as Watson equations [29] ) come from unitarity and crossing symmetry 
and their explicit expressions are 

F n (6i, . . . , $i, ...,#„) = S(6{ — Oi+i) F n (9i, . . . , 6 i+ x,6i, . . . , 9 n ) , (40a) 

n 

F n {9i + 2m, 9 n ) = [J S(9i - 9 X ) F n {9 x , 9 n ) . (40b) 

i=2 
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FIG. 2: Graphical form of the Watson equations. 



Their graphical representations are given in Fig. [2j 

The recursive equations, on the other hand, come from the pole structure of one-particle intermediate 
states. The form factors of integrable theories have, in general, two kinds of simple poles in the strip 
< ImOij < 2ir and except for these singularities they are analytic in the strip. The first kind of poles 
corresponds to kinematical singularities at (9^ = in and their residues give rise to a set of recursive 
equations between the n-particle and the n + 2-particle form factors (see Fig. [3]) 



i Res F, 



n+2\ 



l )=(i-f[s(e-9 i )j 



(41) 



The second kind of poles is instead related to the bound states of the theory. Since there are no bound 
states in the sh-G model, there are no such poles in the form factors of this theory and we do not need 
to write here the corresponding residue equations. 



The general solution of the Watson equations ( 40 1 can be written as 



F n (9 1 



■,0 n )l[F n 



(42) 



where the factors K n {6\, . . . , 6 n ) are completely symmetric and 27ri-periodic functions in all di, and 
-Fmin(#) is an analytic function in < Im6* < ir (without zeros and poles in this strip) which tends to a 
constant value for large values of 9. This function satisfies the equations 



F min (9) = S{9) F min (-9) . 

F min (iTT - 9) = F min (iTT + 9) , 



(43a) 
(43b) 



and its role is to take care of the monodromy properties of the form factors as ruled by the Watson 
equations. The equations ( |43[ ) and the analyticity requirement of -F m i n are able to fix this function up 
to normalization, as we will see explicitly for the sh-G model. The factors K n (9i, . . . ,9 n ) in ( [42] ) must 
contain all the expected kinematical poles. In addition, they must fulfill the proper asymptotic behavior. 
This yields the following final parametrization of the generic n-particle form factor 



F n 



,9 n ) — Fl n Q n (xi, 



on 



F, 



i<3 



(44) 




FIG. 3: Recursive equation of the kinematical poles. 
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where H n is a normalization factor, Xi = e 9i and Q n (xi, . . . , x n ) is a symmetric polynomial. In view 
of Eq. p8| , the polynomial Q n of a scalar operator has the total degree equal to the degree of the 
polynomial Ili^^i + x j) m the denominator, n{n — l)/2. The actual expression of the polynomials Q n 



can be determined by solving the recursive equations ( 41 1 . To this aim, it is convenient to make use of 



(n) 

the basis given by the elementary symmetric polynomials a k of the n variables Xi defined by 



f[(x + Xi ) = x n ~ k a^ (xi,...,x n ), (45) 



t=i fe=i 

or explicitly 



^ Xi t ...x in . (46) 



It is worth remarking that there is a one-to-one correspondence between the infinite set of form factors 
{F n (9\, . . . , 9 n ) , n 6 N} which are the solutions of the functional and recursive equations, and the 
operator content of a relativistic field theory (see [41]). Let us see how this correspondence is realized in 
the sh-G model by solving the form factor equations of this theory. 



B. Form factors of the Sinh— Gordon model 



For the sh-G model the minimal form factor F m i n (6), solution of equations (43), is [30] 



F ,,n V ,„ } J / d< sinh (f a ) sinh (|(1 - a)) , f t§ . . , 
F min (0) - M exp { 4 / - - sinhWcosh(|) sin I — ) ) , (47) 




where 9 = iir — 8 and the normalization constant M — F m i n (iir) is chosen to be 

?dt sinh(|a)8inh[|(l-a)] | _ 1 
^~ P ^ i i sinh(t)cosh(§) r-„™^^ ex P 




In addition to the functional equations (43), F min (9) also satisfies 



F min (iTT + 9)F min (0) = . r *? > • (48) 
sinh 6* + smh(i7ra) 



With the choice 

_ /4sin(7ra) 



#2n = #2 



4sin(7ra) N " 



TV 



the recursive equations for the polynomials Q n entering (44) can be written as 

(-l) n Q n+2 (-x, x, xi, . . . , x n ) = xD n (x; xi, . . . , x n )Q n (xi, ...,x n ), (49) 
where the functions D n are given by 

n k 



D n (x;x 1 ,...,x n )=Y: E (-l) (fc+1 W (n -* )+m 4"M"-L- (50) 



fc=l m=l,3,5,... 

In this formula crl a- re the elementary symmetric polynomials ( 46 1 while 



s sin^ro) 
sm(7ra) 



10 



As shown in |31j , a solution of the recursive equations ( 49 1 is given by the class of symmetric polynomials 

Q„(fc)=detM„(fc), (52) 
where M n (k) is a (n — 1) x (n — 1) matrix with elements 



= (J (») 

°2i-j 



(53) 



/ 4sin(7ra) 



with the exponential fields e kg ^ [33 E2]. With the normalization given by H n (k) = y — 
explicit form of all form factors of these operators is then 



i/2 



The corresponding form factors are the matrix elements of a continuous family of operators identified 

k], the 

(54) 



F n (k) = <O|e fc ^|0 l5 2 



[k] 



l<3 



X 2 | X r. 



So, for instance, the one and two-particle form factors are given by 



M<t>\a\ — 



sin(fc7ra) 
77 ^/sin(7ra) ' 
4 sin 2 (fc7ra) 
M sin(7ra) 



F ■ 



(55) 
(56) 



From now on we will concentrate our attention on the form factors of the even powers of the field </> 
because we will need only these operators for the future computation of the expectation values of the LL 
model (only even operators can have non-zero expectation values). It is useful to express the operator 
content of the theory in terms of a class of particular operators, denoted by : cf> k : , which start creating n 
particles out of the vacuum only when n > k: 



,e„ ) = if n < k. 



(57) 



In their form factors for n = k the polynomial term Q2k(%i ■ ■ ■ , %2k) is equal to the polynomial rii^C 2 ^ + 
Xj) of the denominator and they cancel each other, giving 



Ff--{e 1 ,...,e k ) = 2 k k\ 



Mg 2 sin(7ra) 



mm \ w ij 



(58) 



In view of the recursive equations (41 1, the absence of kinematical poles in ' (6±, . . . ,9k) obviously 
implies the vanishing values (57). To compute the form factors of these operators when n > k, we can 



take advantage of the knowledge of the form factors (54 1 of the exponential operators. Let us denote by 



<fi m the operator whose form factors F™ are obtained by extracting the 0{k m ) term in the expansion of 
F n (k). In view of equations (57) and (58) we have 



n n 

which implies a mixing among the operators : 4 

<P 2 = :cj) 2 : 
M — - At- 



fc-2 



(59) 



;=2,4,... 



At 



fc-2 

E A 

l = 2A. 



(60) 



In Appendix A we discuss how to compute iteratively the coefficients Af. 
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C. LeClair— Mussardo formalism 



At equilibrium the expectation value of a local operator 0{x,t) at temperature T and at finite density n 
is given by 

Trie 'bt O) 

(0) T . n = V W • (61) 
Trfe fc B-r 

For translation invariance at equilibrium (0) T is independent of x and t. If we specify this formula to 
an integrable quantum field theory and we use the basis of multiparticle scattering states, we have 



(0> T ,„ = ^-Efci / ^-^(lI e " T ]^ W.")^ fc), (62) 




where ^T,n = Tr I e k B T ) . As shown in |34| . this expression can be neatly written as 



1 f dfli 

A ■ -n 




( )T,n=2.fc! / -^•■•^Hl 1 - ^ lfe..,W,0)|^,...,^ m , (63) 



where e(0) is the pseudo-energy, solution of the Thermodynamical Bethe Ansatz equation of the model 
of interest, while the connected form factor is defined as [33] 



(6 k , 9 1 \0\9' 1 , . . . , 6' k ) CODn = F (Urn (Q\0\6' x , . . . , 9' k , 9 k - iix + i Vk , ...,6 1 -i7r + i m ) , (64) 

where T in front of the expression means taking its finite part, that is, omitting all the terms of the 
form rji/rjj and 1/rjf where p is a positive integer. In Appendix [B] we give an explicit example for the 
calculation of the connected limit. In this formulation the fi- and T-dependence of the right hand side of 

63) is hidden in e{9) that, for the sh-G model, satisfies the fx- and T-dependent equation (29). Expression 

63 



was checked in various cases [HI go] and was compared with the direct evaluation of the expectation 
value (61) using finite volume regularization [46] . 

Notice that for the sh-G model, in view of the functional relation (48) the connected limit (64) for the 
product of the i 7 ' m i n (%) in the form factors of O simply becomes 

2k k 
n F min(%) — > {Fn»n{nr)) k Y[- 



i t— j- sinh 9ij sinh 9ji 



sinhfi, + sinh(i7ra) sinh#.,i + sinh(i7ra) 

i<j i<j J ' J y 

sinh 2 9 



Af k TT g 13 2 . (65) 

j^. sinh 9ij + sinh (na) 



This means that it is not necessary to employ the explicit form (47) of F mm (9) to calculate the connected 
form factors and, for their actual determination, we only have to take the connected limit of the rest 
of the form factor formula. In particular, the connected form factor F 2 £ ' can be calculated from the 



explicit formula (58) for F^t ' ■ Since it only depends on the rapidities through the F m i n factors, using 
(65) we can immediately write down the connected form factor: 

F^ : -2 2k (2kv( — Vn Shlh2 ^ (66) 

( 2k >-{gi s - m (na)) Wsmh^+sinhVc*)' ' 



IV. THE DOUBLE LIMIT OF THE SINH-GORDON MODEL 



In this section we analyze in detail the mapping between the sh-G and the LL models. We show that it is 
possible to obtain the LL model from the sh-G model by taking the non-relativistic limit simultaneously 
with the limit g — ¥ 0, where g is the sh-G coupling constant. In particular, we show how this mapping is 
realized at the level of the S'-matrix, the Lagrangian densities and the Thermodynamical Bethe Ansatz 
equations. 
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A. Double limit of the two-particle S-matrix 



Let us consider the exact iS-matrix of the sh-G model 

sinh(9 — i sin(a7r) 



S sh . G (9,a) = 



sinh# + i sin(a7r) ' 



(67) 



and let us take its non-relativistic limit accompanied by a simultaneous limit of the coupling constant g 
toward smaller values such that 



oo , g — > 0, g c= fixed 



The resulting expression 



S(9,a) 



Mc 



fg 2 c 



p 

Mc 



coincides with the LL S- matrix (11 1 once we set the sh-G and LL masses equal, M = m, and 



2^,2 



A 



K z c 
~A6 



9 



(69) 



(70) 



Hence the S'-matrices of the two models coincide in this double limit. It is worth noticing that the resulting 
coupling A of the LL model does not need to be small and therefore we shall be able to study the LL 
model at arbitrarily large values of its coupling. To use this correspondence between the two models to 
calculate correlation functions in the LL model, we need to establish the relation between the operators 



of these two theories. For this reason in the next section we show how to perform the limit (68) on the 
fields and the Hamiltonians. 



B. Non-relativistic limit at the Lagrangian level 

Consider the sh-G Lagrangian density 



C 



1 / d(j> 



1 



2 2 



(cosh( ff 0) - 1) 



(71) 



2 \cdt J 2 \dx J g 2 h 2 
To study its non-relativistic limit, it is convenient to write initially the real scalar field in the form |47H50j 



(x,t) = 



I h 2 

2m 



c 2 

il>(z, t) e~* * 4 + ^(x,t)e +t ~Tr 



(72) 



Substituting this expression into the Lagrangian (71) and taking the limit c — > 00, we can discard all the 



oscillating terms, that is, terms containing factors e lnm ° c l ht (with n non- vanishing positive or negative 
integers). These terms, in fact, oscillate very rapidly in this limit and average to zero when integrated over 
any small but finite time interval. In more detail, the relativistic canonical momentum can be written as 



U(x,t) = \4>(x,t) 



I h 2 1 

2mo c 2 



ip(x,t) 



^{x^e- 1 ^ 1 ~^{x,t)e +t ^ t ) +0 



This allows us to express ip and ^ in terms of 4> and LI (up to order O (-7)) 

'1 [m^ 



tp(x,t) 



h V 2 ^ ' ^ y/2m,Q 



IL(x,t) , 



(73) 

(74a) 
(74b) 
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It is easy to show that the commutation relation 

[(/>(x,t),n(x',t)] = ihS(x-x') 

implies the following commutation relation for the non-relativistic operators 

[ip{x, t), ip\x' , t)] — S(x — x) . 



Turning to the Lagrangian density, the kinetic term K of ( 71 ) becomes 

h 2 



K 



dip' dip 



2moc 2 dt dt 2m { 



— W + . - [ip t _ _ + 2moC ^ . 



Expanding the formula (19) for m in the combined limit (70) we obtain 



M 



2 M 2 A 2 



1 



3 /i 2 c 2 



(75) 
(76) 

(77) 
(78) 



From the second term of ( 77 ) we see again that the physical masses in the two models should be equal 
m = M and then the first term is of order 1/c 2 and can be dropped in the limit. 

Let us now turn our attention to the interaction term cosh(<?0) in the Lagrangian (71), which is 
equivalent to an infinite series in terms of even powers of the field (p. Expressing (p in terms of the new 
fields ip and ip\ we can use the binomial formula to expand each power <p 2k in terms of these fields. 
Taking into account that the oscillating terms should be dropped, only the symmetric "middle term" of 
the binomial expansion survives from each power. Collecting the combinatorial factors from the different 
expansions of the powers, we arrive at the following series: 



u(4>) 



2 2 2 2 00 i 



g 2 h 2 



g2 h 2 Z_, (2n)! 



,2-2 00 1 
,2*2 



g 2 h 2 ^ (2n)\ \2m 
The n = l term of the series, 



2n\ 
n J 



g 2n ^ n ip n 



1 TUqC 2 



^ (n!) 2 g 2 h 2 \ 2m 



2„2\ « 



h 2 g 



^ n ip n . (79) 



(80) 



exactly cancels the last term of ( 77 1 . The n — 2 term becomes 

h 2 c 2 g 2 



16 



-^ 2 ip 2 — > Xip^tp 2 



which is just the interaction term in the LL Lagrangian. The rest of the series can be organized as 



E 

n— 3 



2 n {n\f 



E 

n— 3 



->3n-4 



A™ 

(n!) 2 (mc 2 



V> f n ip n + ... 



(81) 



(82) 



where the dots indicate possible higher order terms in 1/c. If we now take the limit ( 68 1 all the terms in 
this series vanish because A is fixed while c — > oo. 

In summary, in the double scaling limit the Lagrangian density ( 71 1 of the sh-G model becomes the 
Lagrangian density ^ of the LL model 



C 



c = — 



2m 



(83) 



So by keeping the coefficient of the ip 4 term fixed, which is the actual constraint enforced by the double 
limit (68), all the higher order terms go to zero and we are left with the non-relativistic LL Hamiltonian. 
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C. The non-relativistic limit of the Sinh-Gordon TBA equations 

To study the non-relativistic limit of the TBA equations of Section |II B| it is convenient to make the 
coordinate change (using from now that m = M) 



Using 



this implies 



and 



For the sh-G model 



p = mc sinh 6 , dp = mc cosh dd . (84) 



d0p«(0) = f = / <WV>. <*■>> 



pM(p) = pW (9(p)) « — p« (— ) (86) 

mc cosh tf(p) mc Vmc/ 

*)«-p(- )• (87) 
mc Vmc/ 



, sin(a7r)\ 2 sin(a7r) cosh(#) . . 

X = -2 arctan 2 =^ y> = . 2 V . ^ L , 88 

sinh(0) J smh (0) + snr {air) 



and in the double limit (68) the kernel <p(0) becomes 



AHmX 

<p(0) — > mc p(p) = mc ^ - . (89) 



Therefore the TBA equations transform into equations ( 13 ) 



2np(p) = l+ f dp'^ip-p')^^'), (90a) 



h 



3G 



— oo 

and p/p^ = 1 + e £ , where 

E(p)= E (£-), (91) 
Vmc/ 

ft = fx — mc 2 . (92) 

Observe that it is correct to take the small p and limit in the integrands even though the integrals are 
extended to arbitrarily large momenta: as a matter of fact, the integrals have a finite support because 
their integrands vanish asymptotically very fast. The expressions for the energies become 



E _ E- Nmc 2 
L ~ L 



F F-Nmc 2 . knT 



L L 2nh 



/ dp^pM(p), (93a) 
J 2m 

dp log (l + e - £ ») , (93b) 



OO 

2 
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which coincide with equations (15). Similarly, the limit of the T = equations (33) is given by 

B 



27rpM (p) = i + J dp' <p(p- p') p« (p') , (94a) 



B 

„2 



£ ~o(p) = -A+^+ / ^^(p-j/)eo(|/), (94b) 



once again in agreement with equations ( 16 ) 



We saw at the end of Section |II B| that the pseudo-energy describes the dressed energy of the excitations 
of the system, which is given by 

E{6) = k B Te{6) + fx (95) 

for the sh-G model and by 

E{6) = k B Te( P ) + fi (96) 
for the LL model. It is worth observing the different behaviors of the excitation energies in the two 



models. As it can be seen from the TBA equation (29), the sh-G energy has a gap M that implies that 
the correlation functions decay exponentially. On the contrary, the LL excitation energy starts as p 2 /2m 



for small momenta (see Eq. (90b)), implying a power-law decay for the correlation functions. Our double 



limit takes care of this difference automatically and thus it will give correct results for the LL model. 

V. LOCAL CORRELATORS FOR THE LIEB— LINIGER MODEL 

In this section we calculate LL one-point correlation functions at fixed particle density n and temperature 



T by applying the formulas of Sections O and III in connection with the double limit presented in 



Section [IVB| The fields are taken at the same position and time: since our system is taken at equilibrium 
and translationally invariant, their correlators are obviously space and time independent. We focus our 
attention on the local fc-particle correlation functions gk defined as 

(^Vl^'fttv) , (97) 

where 7 and r are given in ^ and ^ . These local correlators play an important role in experiments with 
ultracold bosons since the pair correlations are responsible for the rates of inelastic collisional processes. 
Furthermore, the low-temperature recombination rate for a Bose gas is proportional to the local three- 
body correlation function |51) and measurements of the three-body recombination rate can be used to 
determine the local correlations and as a tool for distinguishing condensed and non-condensed phases 
[52] . For 32(7, t) and #3(7, r = 0) exact results are available [20] [25l [26], whereas for the others the 
asymptotic behavior of the correlators in the regimes of small and large coupling or temperature was 
computed in [501 HI] ■ We will present a comparison with the exact and the approximate results present in 
the literature, showing the improvement that our method brings in the computation of these correlators. 
Let us start with the correspondence between the sh-G and the LL operators 

f^-\ k f 2k ](^"^), (98) 



2m J \k 

which can be established along the limit procedure described in Sect ion |IVB| To compute the expectation 
values of these operators at finite density n of the Bose gas and at finite temperature, we need to employ 



the results of Section III so that 



where f(9) = 1/(1 + e e ^) are the filling fractions and the notation "lim" denotes the double limit (68) 



Note that in Eq. ( 99 1 the terms with I < k are zero and therefore the first non-zero term in the series is 
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a fc-fold integral. Now, similarly to what happens for the TBA equations, the filling fractions effectively 
cut off the integrands at large values of the rapidities, so we can exchange the order of the limit and the 
integrals, arriving at a fully non-relativistic formula 



k ip k 



-k 



1 [ dpi , f dpi 



l=k 



iPlJc 



(100) 



Here f(p) — 1/(1 + e E ^) where e(p) is the solution of the non-relativistic TBA equations ( 90|85 ) and 



Fit : (fe»c 



= lim 



inc. 



mc 



(101) 



are the double limit of the connected form factors. We go through the steps of the calculation of a specific 
form factor and we list the explicit expressions of the first few of them in Appendix [B] 



A first check of the validity of formula (100) is provided by the correlator (tjj^(x,i)ijj(x,t)). With the 
explicit connected form factors of : (f> 2 : the first terms read 



t,/,\ - 



tr / / $/(^)/^)/fe)^(pi 2 )^2 3 ) + ..., (102) 



where we use the notation pij = pi — pj and (p{p) is the scattering phase shift in the Lieb-Liniger model 
(14 1. The pattern in (102 1 persists for the further multiple integrals and one easily recognizes that the 



right hand side of ( 102 ) is nothing else but the recursive expansion of 



J dpp (r) (p), 



(103) 



where (p) is the iterative solution of the integral equation ( 90a I 



rW r) (p) 



2irh 



dp' 
2tt 



f>(p-p')p (r) (p') 



(104) 



In this way we successfully recover the identity (tp^tp) — n - While the result may appear obvious, it is 
worth stressing that it was obtained by taking the double scaling limit of the sh-G form factor expansion 
and employing the LeClair-Mussardo formalism, so it provides an important check of the method. 



A. Correlators at T = 



At zero temperature, similarly to the TBA equations, the formula ( 100 ) for the expectation values admits 
a simpler expression 



1 / ,_9 \ —k oo 

E 



k J \2m) 



1 f dpi 



l=k 



11 J 2tt 

-B 



dpi ^:0 2fc :/ , 
ITT 



(105) 



-B 



where B is the Fermi momentum determined by the set of TBA equations ( 94 ) together with 

B 

ii = l p {r) {p)dp. (106) 
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The equations become more transparent by introducing the dimensionless quantities 



, v(k) = hp {r) (Bk), 



2m X hwy 



h B 



B 



(107) 



where in the last expression we used the definition of the LL parameter 7 of Eq. ([2|. In terms of these 
new variables, equations ( 106|94a ) become 



1 

1 = 1/ "{k)dk, 



u{k) 



2n 



20 



dk' 



2tt (k-k') 2 +0 : 



v{k') , 



(108a) 
(108b) 



while the (non-relativistic) ground state energy is given by 



E 
L 



—B 



2m 2m 



,3 (1 

P 



dkv{k)k 2 



2m 



n 3 e(7) . 



(109) 



The strong coupling expansion of (108a I is obtained by plugging into it the iterative solution of (108b I 
and expanding the integrals in 0~ x 



+ ^/r 2 (i - V 2 + V 4 + ..-) + -V 3 (i - 2 ^ 2 + ■•■) + 

IT^ 3 3 7T 4 

This provides a series expansion relation between and 7 



7 = ^-2+ „ , 



4 

3~0 2 



P=- 7 + 2- 



4^ 



which is equivalent to a Fermi momentum 



2 4 

5 = firi7r [1 ^ — 

7 7 



(110) 



(111) 



(112) 



Using the formulas above, we can now derive the leading order behavior in 7 _1 of the general correlator 
gk(l)- It is easy to see that the leading order comes from the first non-zero integral in the series (1051 
where the integrand is the double limit of the connected form factor (66). Taking the double limit, for 
the first term of (ijr k ip k ) we get 



2m 



B 


B 


f dpi 




J 2tt 


■J 


— B 


-B 







^2 2fe (2fc)! ( — 
2tt y ' \8m 



k k 

n 



h 2 P 2 



: fc! 
(27r) fe 



dfci. 



h 2 p 2 - + 4m 2 A 2 



(113) 



For the leading order behavior of these quantitie s, w e nee d to keep only the leading order term 
Qg-n(n-i)) of the i nte g ra nd and substitute from ( (ll2| and ( [lTT| £ = firm and /3 = 7/77. The final 
result is 



_ fc! 
.9fe - 7^ 



fc(fe-i) 



(114) 



18 



0.1. 



0.01 



0.001 : 



1 10 20 

y 

FIG. 4: Deviations 1 — gi from the exact result (gi = 1) at T = using form factors up to n = 
respectively with green dot-dashed, blue dashed and red dotted lines. 



6 and 8 particles, 



where 



1 1 k 

I n = j dk x ... fdk k f[, 

-1 -1 1 <J 



(115) 



So, in this way we recover the expression obtained in |20j by using completely different methods. Let us 
now discuss in more detail the results for g%, gi and g%. 



1. Expectation value g± 



We already demonstrated that our series expansion sums up to the exact value gi = 1, here we show how 
convergent the series is. The actual computation consists of the following steps. First we solve numerically 
the integral equation (108b I, second we obtain the /3(7) function from (108a), and finally, we integrate 



numerically the dimcnsionless forms of the integrals on the right hand side of (105 1. 

Carrying out these steps for gi = (ip^^/n, we obtain the plot shown in Fig' [4] for g\ as a function of 
7. The series ( 105 ) is nicely saturated by the first few terms for sufficiently large values of 7 (one should 



keep in mind that 7 = is a singular point of the LL model and therefore one cannot expect a priori 
any fast convergence nearby). In fact, this plot shows that the exact value gi = 1 is rapidly approached 
by just the first terms of ( 105 1. It is clear that including more terms in the series, (i.e., employing higher 
particle form factors) extends the fast convergence toward smaller values of 7. Notice, however, that the 
convergence of the series is always remarkably fast for all 7 > 1.5, where the exact value is obtained 
within a 5% accuracy just using its first four terms. 



2. Expectation value 32(7) 

Let us continue our discussion with the calculation of the correlation function (ip'ip'ifjip) = n 2 gi. This 
correlator is the expectation value of the interaction term in the Lieb-Liniger Hamiltonian Q, thus it 
can be exactly determined [30] US] via the Hellmann-Feynman theorem [53] : 

, ,+ ,+ v 1 /dH\ d f E\ de( 7 ) , . 

^ t ^ = x(dA) = dAU 52 = ^T' (U6) 



where e(^f) was defined in (109 1. Our result can be compared with this expression, providing a good 



possibility to check again the correctness of the approach. We have 

— B -B -B -B -B 
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FIG. 5: Plot of <72 as a function of 7 at T = using form factors up to n = 4, 6 and 8 particles, respectively 
with green dot-dashed, blue dashed and red dotted lines. The exact value is given by the solid li ne w hereas the 
dot-dot-dashed line below, indicated by the arrow, corresponds to the strong coupling expansion (1191. 



In terms of the dimensionless quantities the expansion in (3 up to the four-integral term gives 



(3 2 \37T 2 ' 



24 



3tt 3 



o 

3-2 



332 
105 y 



3tt 4P \ 15 P 



Substituting the relation (111) we obtain the strong coupling expansion 



.92 



4tt2 

3^ 



6 

7 



(24- 



1 



(118) 



(119) 



The leading order behavior agrees with ( 114[ ) but it is worth noticing that we also obtained subleading 



terms in 7 . We notice that this result can be obtained using the Hellmann-Feynman theorem and the 
expansion of the ground-state energy given in |54] . 

The plot of 32(7) _ obtained by numerical integration and using the integral equations (108 1 for j3{^j) 



- is drawn in Fig. [5j As for the previous example, we see that by increasing the number of form factors 
employed in our series, our result rapidly converges to the exact value. The discrepancy between the exact 
value and the one obtained with four integrals is less than 3% for 7 > 2. Expression ( 119 1 is also plotted 
in Fig. [5] to show that the determination of 52 (at finite 7) obtained from the first terms of Eq. ( 105 1 
is much closer to the exact result. This is because we solve the TBA equations (108 1 with an arbitrary 
precision and every term of our series contains infinitely many powers of 7. 



3. Expectation value 93(7) 

As a final example let us discuss g 3 , a quantity known exactly up to now only at T = [55] . An analysis 
similar to the previous cases reveals that 



3 r 



16 



/3 3 L15tt 



■0- 



144 „ , 



32 



/r 7 (1 + ...) + 



(120) 



Trading j3 for 7 using (111) we arrive at 



33 



16 7T 6 

15 76 



16 

7 



+ G(7~ 8 ). 



(121) 



Here the leading order term is the asymptotic result (114), but as in the previous example, we also 
obtained the next order in the large 7 expansion. 
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0.0001 



FIG. 6: 33 at T = with form factors up to n = 6 and 8 particles with blue dashed and red dotted lines, 
respectively. The exact value is given by the solid line whereas the purple dot-dot-dashed line above corresponds 
to the leading order expression (|114[). 



The logarithmic plot of g% using the form factor expansion up to n — 6 and 8 particles (one or two terms 
from the series) is shown in Fig. [6] together with the exact result of [55]. As in the previous examples, 
this plot shows a nice convergent pattern toward the exact value. The leading order (114 1 in the large 
7 expansion is also plotted in Fig. [6] to show that in this domain of 7 this result largely differs from 
the exact value. The subleading term, given in (121) provides an improvement for larger 7, however, the 



result is still quite far from the one obtained with our method. 



B. Correlators at finite temperature 



To obtain the expectation value at finite values of the temperature T we have to employ the formula 



100 1 which contains non-trivial filling fractions and we need to solve the whole set of TBA equations 
85|90 1 . It proves to be useful to introduce now a different set of dimensionless quantities 



which satisfy 



?=-?-, a =y4^: 9(1) = - p( nh i q) 1 ( 122 ) 



M = -a + ^ - / ^ , * 1 log (l + e-W>) , (123a) 
r J 2tt (q — q') z + 1 \ / 

—00 

00 

g(l) = 77- h [ dq' -. %— - g{q ' ] , ,. , (123b) 

y j 2ira J (q - q') 2 + 1 1 + e £ (? ') ' v ; 



and 



J_ = f 9(g) 



dq. (123c) 



The physical parameters of the problem are A, T and n, but only the dimensionless combinations 7 and 
t enter the results. The chemical potential (or the dimensionless fugacity-like parameter a) gets fixed by 



the constraint given by the last equation. Rescaling the density with a in ( 122 1 makes it possible to find 



the self-consistent solution by an iteration of the system of equations. Once we have e(q) for a given 7 



and r we can substitute it into the form factor expansion (100). 

Before discussing the numerical results, let us derive an asymptotic formula in the regime where 1 <C 
t <C 7 2 . The second inequality implies that the q 2 term in the equation for e(q) becomes quickly quite 
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FIG. 7: Deviations 1 — 51 from the exact result (g\ = 1) as a function of the scaled temperature r for a fixed 
value of 7 = 7. Inset: 1 — g± vs 7 at r = 1. In both figures form factors are used up to n = ! (green dot-dashed), 
6 (blue dashed) and 8 (red dotted) particles. 



large while the first one implies that a is a large negative number. This means that, even for q close to 
zero, the convolution term is small and the leading order e(q) is given by 



e{q) = -a - 



2 9 

97 



(124) 



so we can make an expansion in the small parameter exp(— a + q 2 ^ 2 /t) [55] , At the leading order 
g{q) = l/(27ra) and this implies 



dq a _<W 1 
— e T = — 
2tt 7 



(125) 



so we arrive at the 7-independent result 



(126) 



We see that for r > 1 a is indeed a large negative quantity. 

Due to the condition 7 3> 1, we can again restrict ourselves to the large 7 limit of the first non-zero 
term in the series (100). Substituting the Boltzmann filling fraction f(ff) = e~ € given above, we arrive at 



where 



Jk ■ 



(127) 



Jk 



k\ 



dxi . . . dxk e 



IK 



2 fe (fc-i)/2 



(128) 



with Bk+i = (k + l)r(fc + 2)Sfc, B\ = 1. This is exactly the result found in [20]- However, as one can check 
numerically, the filling fraction comes close to a Boltzmann distribution only for such extreme parameter 
values as 7 ~ 1000 and r - 10000. In Figs. |] and [9] one can see how large the difference is between this 
leading order approximation and our result for r = 10. The improvement achieved in the determination 
of this quantity by the method proposed in this paper may have an important experimental relevance. 
Let us turn now to the numerical results obtained by exactly solving the TBA equations (123), sub- 



stituting e in the formula ( 100 ) and numerically integrating the first terms in the series. We consider 
separately the computation of gi, <? 2 and g^. 
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FIG. 8: J2 at t = 1, 10 using form factors up to n = 6 and 8 particles with green dot-dashed, blue dashed and 
red dotted lines, respecti vely. The solid lines show the exact result, while the purple dot-dot-dashed line is the 
leading order expression (127 1. 



1. Expectation value g± 

To test the reliability of our expansion at finite temperature we computed the deviations from the exact 
result for the trivial expectation value gx = 1. The results are plotted in Fig.[7J showing that the precision 
does not decrease with increasing temperature: even using only three terms of the expansion (i.e. summing 
up to n = 6 particles) the error is < 1% in the range of temperature between r = and r = 15. In the 
inset of Fig. [7] we plot the deviation from the exact result as a function of the LL parameter 7 at a fixed 
temperature. 



2. Expectation value (72(7,1") 



For T > the Hellmann-Feynman theorem gives 



d If 

dA i L 



(129) 



where the free energy can be calculated from the TBA approach ( 93b ) . In dimensionless variables 



92(l, t) 



d7 



dq 
2tt 



log(l + e 



(130) 



We derive now a simpler expression for this. The trick is to substitute for 1/27T under the integral the 
rest of Eq. ( 123b ), then using the associativity of the convolution by an even function and finally use the 



derivative of Eq. (123a) with respect to 7. Many terms drop out and we are left with 



.92 = 2 7 2 



dq 



Qjg(g) 
1 + e^te) 



dq 
2tt 



log(l 



(131) 



The advantage of this expression is that it is enough to solve the TBA equations for the value of 7 we are 
interested in instead of evaluating the free energy for several 7 and then differentiating it numerically. 

Our evaluations of 52 at r = 1 and r = 10 based on the form factor expansion are shown in Fig. [8] 
together with the exact result (131) and the leading order result (127). The convergence of our series 



is basically as good as it was for T 
7 < 10. 



and it is clear that the asymptotic formula fails, especially for 



23 



S 3 

1 
0.1 

0.01 
0.001 
0.0001 



T=10 



%=1 



10 



20 



FIG. 9: gs, vs 7 at r = 1 and r = 10. The blue dashed and the red dotted lines refer to n — 6 and 8 particles, 
respectively; the purple dot-dot-dashed line shows the asymptotic result (127 1. 



3. Expectation value 513(7 , t) 



In the case of 52 and of 33 at T = we checked our results using exact formulae. We learned that on the 
one hand our result for 173 reaches the same accuracy as 172 for slightly higher values of 7, on the other 
hand going to finite temperature does not spoil the precision of our results. Thus we are confident about 
the reliability of our results for g% at T > 0, at least for not too small values of 7. We emphasize that 
neither exact results nor approximations of precision comparable to ours exist in this case, which renders 
our evaluation a new result. 

Fig. [9] sho ws g% as a function of 7 at fixed temperatures t = 1 and r = 10. In the figure the asymptotic 
result (127 1, va lid for large temperature and coupling, is also plotted: one can see that even for r = 10 
the result (1271 does not give the exact asymptotic behavior (which is only reached for very large values 
of the scaled temperature r). Fig. 10 shows instead g% as a function of r at fixed values of 7 = 7 and 
7 = 15. The asymptotic formula ( 127) is different from our result by a factor of ~ 10 at 7 = 7 and r = 10. 



VI. CONCLUSIONS 



In this paper we have shown that the Lieb-Liniger model, describing one-dimensional interacting bosons, 
can be obtained as a non-relativistic limit of an integrable relativistic field theory, the Sinh-Gordon 
model. In this limit, the S'-matrix, the Lagrangian and the Thermodynamical Bethe Ansatz equations of 
the sh-G model reduce to those of the LL model. We have also shown that the pseudo-energies of the 
sh-G TBA (which are actually the energies of the excitations above the vacuum) become massless modes 
in the non-relativistic limit, in agreement with the hydro-dynamical description of the LL model given 
by bosonization. 

The mapping between the two models proved to provide an efficient method to compute expectation 
values in the LL model by using the form factor expansion of the expectation values of the relativistic 
counterpart. The main advantage of using the form factors of the relativistic integrable sh-G model is 
that, as for any relativistic quantum field theory, its form factors obey a set of stringent constraints 
that permits the determination of their exact expressions. Moreover, the quantum integrability of the 
relativistic model allowed us to employ the rich collection of results valid for these systems, like the TBA 
and the LeClair-Mussardo formalism. 

Using these two formalisms of the relativistic sh-G model (form factors and TBA), we computed the 
expectation values of the LL model. The method works equally well at T = and T / where the 
series expansion presents a remarkable convergence behavior for finite values of the LL parameter 7. The 
computation of one-point correlators was presented in detail, as well as the comparison with the known 
results available in the literature. In particular, we have determined the expectation value 33(7, r) at finite 
temperature, for which there had been only asymptotic analytic results in the literature. This quantity 
is related to the recombination rate of the atomic gas and thus to the lifetime of the experiments. 

It would be interesting to analyze the possibility of extending our method to other non-relativistic 
strongly correlated systems, identifying their relativistic counterparts. An equally important direction 
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FIG. 10: 33 vs the scaled temperature r for 7 = 7 and 7 = 15. The blue dashed and the red dotted lines refer to 
n — 6 and 8 particles, respectively; the purple dot-dot-dashed lines show the asymptotic result (1271. 



would be to see how the methods presented here can be used to compute space- and time-dependent 
two-point functions. 
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Appendix A: Determination of the coefficients A\ 



To find the coefficients A k at level k in Eq. ( 59 ) we can proceed in the following way. If we already know 



the form factors F,' ' for j < k, we can extract the A k, s one by one from Fh 



F k 



F 



Ft 
F 2 ' 



A k = 



Pf - A k F 4 q 



(Ala) 
(Alb) 



Af 



fc - 2 Ak jp-> 

3=2,4,- A j r l 



Once these coefficients are known, all the form factors of : <j> k : , including the ones with n > k 
can be obtained from F„. These will be needed at higher levels. For example 



F 



10 



F 6 



46 F -4>' 



F 6 — A 6 

r 10 ^4 



A 2 r 10 



A 6 F 2 

A 2 r W 



As an example we give the rules for the "operator mixing" ( 60 ) explicitly at the first four even 

4.2 . 



A2 - 



A4 - 



2^,2 



ir a 



2„,2 



20 



7T a 



56 



16 



336 



4 4 
7T a 



64 



(Ale) 
particles 

(A2) 
levels: 
(A3a) 
(A3b) 

(A3c) 

(A3d) 
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Appendix B: Explicit formulas for the connected form factors 



In this Appendix we explicitly go through the steps required for the calculation of the limit of the form 

~ ,2k , . . 

factors F 2 [ ' ({pi}) CO nn that enter into the formula (100), and then we list the first few of them. 

We consider explicitly only the simplest non-trivial case, ' (pi,P2)conn- We start from the four- 
particle form factor of the exponential operator (54 1: 



F 4 (k) = (O|e fcff *|0i,0 2 ,0 3 ,0 4 > 

[k] 



4sin(7ra) 



[kfvw* -[k- l][k][k + 1}(<J 2 3 - a 2 a 4 )) f[ Fmin( ^ ^ 



(Bl) 



where the elementary symmetric polynomials ( 46 ) are to be understood as <7k = cr[ 4 ^ and we recall that 



Xi = exp(6 l i) and [k] = sin(fc7ra)/ sin(7ro!). The form factor of :<fr 2 : is given by the 0(k 2 ) term in the 
Taylor-expansion in k (for : 4> 2 : there is no mixing ( 60 1 ) : 



F; 



'1, C2, 03, u 4 



32n 2 a 2 2 2 jj 

y i<3 



e 6i _|_ e 8j 



(B2) 



To perform the connected limit ( 64 1 we first recall that the product of the minimal form factors in this 
limit is given by (65 1. For the rest of the formula we write x% — ~X2 + M?2; x 4 = —x\ +irji and we expand 
the polynomials in the numerator and denominator to obtain the finite part: 

x\{x\ + x 2 2 ) rft + 4x\x\ mm + x\(x\ + x\) r$ + C(t7 3 )" 



T 



lim 

?7i->0,r)2- 



= T 



\x\x\ 



01 - 2x 2 x 2 + x2*) mrm + 0(t] 3 
x%(xI + x%)tii x\{xI + xI)t}2 



_(xj-xl) 2 (X 2 -X 2 ) 2 T]2 

Collecting all the terms we obtain 

^4 [t>l, f2jcom 



{x\-X 2 ) 2 T]l 



O 



(- 






mj. 



\x\x 



2^2 
2 



1 



(x 2 -x 2 ) 2 



sinh 2 ( 



Af 2 sinh 2 6»i2 



M 2 g 2 sinh 2 9 12 sinh 2 6 12 + sinh 2 (7ra) 



327r 2 a 2 
~7 



sinh 2 812 + sinh 2 (7ra) 



(B3) 



(B4) 



We note that another way to arrive at this result is to calculate first the connected form factors of 
exp(kg<f>) and then to extract the 0{k 2 ) term to obtain the connected form factors of : (fx 2 : . 

Finally, after substituting 0i — > pi /mc and using the definition of a ( 21 ) , we can perform the double limit 
( 68|70 ). We list below the result together with the first few non-relativistic form factors obtained in this 
way. We use the notation = pi —pj, while denotes a sum over permutations of {i, j, k} = {1, 2, 3}. 



~ ,2. 

F 6 ' (Pl,P2,P3) 



1 

h' 

8mA _ 2 . 
4m 2 A 2 + h 2 p 2 2 ~ h^ Pl2 >> 

32WA 2 (12m 2 A 2 + h 2 (p 2 2 + p\ 3 + p 2 23 )) 
(4m 2 A 2 + h 2 p 2 2 )(4m 2 \ 2 + h 2 p 2 13 )(4m 2 X 2 + h 2 p 2 3 ) 



(B5a) 
(B5b) 

= \ H ^(Pij)^(Pjk) ■ ( B5c ) 



^4 1 (Pl,P2)conn = — TT ^{P\2)Px2 > ( B6a ) 



F 6 {p\ , P2;P:i)cotm - 8/imA ( . ; 



(gig +P 2 13 +P 2 23 )(8m 2 \ 2 + h 2 {p 2 2 +p 2 13 +p 2 23 )) = 
(4m 2 A 2 + h 2 p\ 2 )(4m 2 \ 2 + h 2 p 2 l3 )(Am 2 \ 2 + h 2 p 2 23 ) 

-rr fiiPij^iPi^Pik ■ ( B6b ) 
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2 2 2 

g 

= 16m 3 A 3 ^(Pl2MP23)<^(Pl3)Pl2 P23 Pl3 • ( B7 ) 
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